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Abstract 

We discuss the possibilities for the direct detection of neutralino dark matter via elas- 
tic scattering in variants of the minimal supersymmetric extension of the Standard Model 
(MSSM) with non-universal supersymmetry-breaking contributions to the Higgs masses, 
which may be either equal (NUHMl) or independent (NUHM2). We compare the ranges 
found in the NUHMl and NUHM2 with that found in the MSSM with universal supersymme- 
try-breaking contributions to all scalar masses, the CMSSM. We find that both the NUHMl 
and NUHM2 offer the possibility of larger spin-independent dark matter scattering cross 
sections than in the CMSSM for larger neutralino masses, since they allow the density of 
heavier neutralinos with large Higgsino components to fall within the allowed range by as- 
trophysics. The NUHMl and NUHM2 also offer more possibilities than the CMSSM for 
small cross sections for lower neutralino masses, since they may be suppressed by scalar and 
pseudoscalar Higgs masses that are larger than in the CMSSM. 
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1 Introduction 



Supersymmetry is one example of a theory that suggests the possibiUty of observing new 
particles at high-energy particle colliders such as the LHC. Many super symmetric models 
predict that the lightest supersymmetric particle (LSP) should be stable and hence present in 
the Universe as a relic from the Big Bang [1] . As such, the LSP may provide all or some of the 
cold dark matter postulated by astrophysicists and cosmologists. In many supersymmetric 
models, the LSP is the lightest neutralino, x, though there are other possibilities such as the 
gravitino. 

The most convincing vindication of any model of dark matter would be the direct detec- 
tion of astrophysical dark matter particles via their scattering on nuclei in low-background 
experiments underground [2]. The model could then be further verified by comparing the 
direct dark matter detection rate with theoretical calculations based on measurements of 
the model parameters at colliders. If the LSP is indeed the lightest neutralino x, many 
estimates suggest that supersymmetric dark matter could be detected directly in present or 
forthcoming experiments [3-14]. On the other hand, if the LSP were the gravitino, the cross 
section for its scattering would be so small as to be undetectable. 

Within the neutralino LSP scenario, the absence to date of supersymmetric signals in ac- 
celerator experiments imposes constraints on the possible direct detection rates. In general, 
these depend on the specific supersymmetric model, and specifically on the way in which 
supersymmetry is broken. We consider here the minimal supersymmetric extension of the 
Standard Model (MSSM) [15], and focus on three possibilities for the pattern of supersym- 
metry breaking. 

In the simplest model the super symmetry-breaking scalar masses mo are constrained to 
be universal at some input GUT scale, a scenario called the CMSSM [16, 17]. The absence 
of fiavour-changing interactions beyond the Standard Model motivates the assumption that 
squarks or sleptons with the same quantum numbers should be (very nearly) degenerate 
[18], and universality between (some) squarks with different different quantum numbers and 
sleptons is suggested by some GUT models. These also suggest that the SU(3), SU(2) and 
U(l) gaugino masses mi/2 should be degenerate. The direct detection of neutralino dark 
matter within this CMSSM framework has been studied in many papers [4,5,7,9-14]. 

However, there are no good reasons to expect universality for the supersymmetry-breaking 
contributions to the scalar masses of the MSSM Higgs multiplets. Accordingly, our main fo- 
cus here is on variants of the MSSM with non-universal Higgs masses [19], generically called 
NUHM models. The supersymmetry-breaking contributions to the two Higgs masses might 
either be identical, a scenario called the NUHMl [20,21], or might be different, a scenario we 
call the NUHM2 [22,23]. The direct detection of neutralino dark matter within the NUHM2 
has also been studied in several papers [6,13,20], but direct detection within the NUHMl 
has been less studied [20]. 

Since the amplitude for spin-independent dark matter scattering on a heavy nucleus 
receives important contributions from the exchanges of Higgs bosons, and since these are less 
constrained in NUHM models than in the CMSSM, a priori one would expect that the direct 
dark matter scattering cross section should exhibit more variability than in the CMSSM. In 
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particular, whereas the present experimental upper limits on dark matter scattering barely 
touch the range expected within the CMSSM, ongoing direct dark matter searches might 
already be sensitive to significant samples of NUHM models. Conversely, whereas the planned 
direct dark matter searches should be able to explore much of the CMSSM parameter space, 
it might be more likely that NUHM models could escape detection. 

In this paper, we first review in Section 2 the relationships between the CMSSM, NUHMl 
and NUHM2 models, as well as basic formulae underlying the calculation of the spin- 
independent and -dependent dark matter scattering cross sections. We then review briefiy 
in Section 3 the expectations for direct dark matter detection within the CMSSM. Subse- 
quently, we explore in Sections 4 and 5 the possible scattering cross sections in the NUHMl 
and the NUHM2, respectively. These explorations each proceed in two steps: surveys of 
the variation in the direct dark matter cross section across some representative planes for 
fixed values of the other parameters, followed by a presentation and discussion of the overall 
range of possible values throughout the corresponding parameter space. We concentrate on 
spin- independent dark matter scattering because this seems closer to the present and likely 
experimental sensitivity than does spin-dependent scattering. Our conclusions are summa- 
rized in Section 6, where we comment on the possibilities for distinguishing experimentally 
between the CMSSM, NUHMl and NUHM2. 



2 Theoretical Framework 

2.1 From the CMSSM to the NUHMl and the NUHM2 

The CMSSM has four continuous parameters and a sign that determine the weak-scale ob- 
servables: the supposedly universal scalar mass mo, gaugino mass mi/2 and universal trilinear 
coupling Aq] the ratio of the Higgs vacuum expectation values, tan/?, and the sign of the 
Higgs mixing parameter fi. The parameters mo, mi/2 and Ao and the sign of fi are assumed 
to be specified at the GUT scale. A useful starting-point for analyzing the CMSSM is the 
(mi/2, mo) plane for some fixed value of Aq, tan /3 and sign of /i. The input supersymmetry- 
breaking contributions to the effective Higgs masses-squared, mf 2 for the Higgs fields coupled 
to the down- and up-type quarks respectively, are renormalized differently below the GUT 
scale, because they are sensitive to different Yukawa couplings, specifically the t-quark cou- 
pling in the case of H2, and this renormalization permits electroweak symmetry breaking 
at low energies. In the CMSSM, and m^ are calculated from the electroweak vacuum 
conditions: 

mliQ) = mjiQ) + mHQ) + 2fi\Q) + A^iQ) (1) 

and 

.2 /I _ o\ I a(1) 



mt 



ml tan^ f3 + |m|(l - tan^ (3) + A} 



^ tan^/^-l + Af ' 



where and A[/'^^ are loop corrections [24-26], Q = (m^^^m^^)^/^, and all quantities in 
(|2| are defined at the electroweak scale, mz- Unless otherwise noted, uia = mA{Q) and 
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/X = jj(mz)- The values of the parameters in Q and (|2| are related via 

ml{Q) =m\ + ci, 

ml{g)=ml + C2, (3) 

where ci, C2 and are well-known radiative corrections [24,27,28]. 

The restriction mi{MGUT) = rniiMcuT) is retained in the NUHMl, but these may differ 
from the otherwise universal scalar mass, itiq. The common GUT-scale value of the Higgs 
masses-squared therefore introduces an additional parameter in the NUHMl. This can be 
mapped into an additional free parameter in the effective low-energy theory, that may be 
taken to be either or niA- As in the CMSSM, one uses ([T]) - (|3| to calculate the relationship 
between m\ and m\ at the weak scale that is required to respect the electroweak boundary 
conditions, allowing for the weakened universality condition at Mqut- If ''^a is taken to be 
the additional free low-energy parameter, then at we have [21] 

mliiaxi^ /? + 1 + Ajf)) = mliisn^ p + mKtan^ - 1) - 

+ {m\ - (A^(g) + ci + C2 + 2c^)) (tan^ 13-1 + (4) 

Alternatively, if /i is taken as the free parameter, then at mz we have 

ml = m^tan2/5 + /i2(tan2/5 - 1 + A^f)) + ^mKtan^ /5 - 1) - A^. (5) 

In each case, the NUHMl boundary condition at Mqut is ml = m\. Clearly, for some specific 
input values of and mA, one finds mKMcur) = mKMcur) = ^o, thereby recovering the 
CMSSM. The subject of Section 4 is the deviation from the CMSSM prediction for dark 
matter scattering that may be found in the NUHMl. We discuss first planar subspaces of 
the NUHMl in which m^ is taken as the additional free low-energy parameter, and then 
consider planes with /i as a free parameter. 

In the NUHM2, the soft supersymmetry-breaking contributions to both the Higgs scalar 
masses mi^2 are regarded as free parameters. Correspondingly, these may be traded for free 
values of both fj, and at low energies. In this case, we can write [22] 

ml{l + tan^ /3) = m\{Q) t&n^ /3 - fi^itan^ /3 + 1 - A^^^) - {ci + C2 + 2Cf,)tan^ /3 

- Aa(Q) tan^ p - ^ml^{l - tan^ p) - A« (6) 

and 

ml(l + tan^P) = m\{Q) - fi\tan^ p + 1 + A^^^) - {ci + C2 + 2c^) 

-Aa{Q) + ^mL(l -tan^/^) + A«. (7) 

Each point in a CMSSM (mi/2,mo) plane can be 'blown up' into a (/x, m^) plane. The 
NUHMl subspace may be represented as a line in such a (//, m^) plane, and the CMSSM as 
one or two points on this line. 
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2.2 Dark Matter Scattering 

The neutralino LSP x is the lowest-mass eigenstate combination of the Bino B, Wino W and 
neutral Higgsinos Hi 2- Its mass matrix is diagonalized by a matrix Z: diag(m^^ ^) = 
Z*NZ~^ , and the composition of the lightest neutralino may be written as 

X = Z^iB + Z^2W + Z^si^i + ^x4^2. (8) 

The amplitude for x scattering on quarks depends on the squark mass-squared matrix, which 
is diagonalized by a matrix rj: diag{m\, m"^) = rjM'^ri'^, which can be parameterized for each 
flavour / by an angle The diagonalization matrix can be written as 

cos Of sin6'je*'^-f 
— sin 6'/e~*'^-'' cos6'j 

The magnitudes of /i and the pseudoscalar Higgs mass uia are calculated using the appro- 
priate electroweak vacuum conditions in the CMSSM, NUHMl and NUHM2, as discussed 
in the previous Subsection. 

The effective four-fermion Lagrangian relevant for relic dark matter scattering is [29]: 

^ = a2tXYl^XQn^Ll^(lt + asiXXQiQi- (10) 

This Lagrangian is to be summed over the quark generations, and the subscript i labels 
up-type quarks (z = 1) and down-type quarks {i = 2). The coefficients are given by: 



) - 




Vl2 \ 




V V21 


V22 J 



(9) 



r Ury ,2 ,^ ,21 T- 



(11) 



and 



4m|cos2^H/ LI X3I I X4I J 2 
1 



where 



+ D^Q (-^ - Re [6u {gZ^2 - g'Z^,)] |, (12) 



x5-i 



2mwBi 



V. ^ ^;i(f^7%i+^?T3.Z,2)+r^;2^|^, (13) 



^We ignore here all possible CP-violating phases. 
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where yi,Tsi denote hypercharge and isospin, and 



6ii = Z^3{Z^4), 62i = Z^4{-Z^3), (14) 

i?j = sin/9 (cos/3), Cj= sin a(cos a), = cosa (— sina) (15) 

for up (down) type quarks. We denote by < niHi the masses of the two neutral scalar 
Higgs bosons, and a denotes the neutral Higgs boson mixing angle. 

In the NUHMl and NUHM2, the greater freedom in the choice(s) of the input values of 
ml and m2 translates into greater freedom for ruHz and rriHi than in the CMSSM and hence, 
a priori, more variation in the dark matter scattering amplitude. 

The elastic cross section for x scattering on a nucleus can be decomposed into a scalar 
(spin- independent) part obtained from a^i (12), and a spin-dependent part obtained from 

(11). Each of these can be written in terms of the cross sections for elastic scattering off 
individual nucleons, as we now review. 

The scalar, or spin-independent (SI), part of the cross section can be written as[^ 

as, = ^[Zf, + iA-Z)U]\ (16) 

71 

where rrir is the x-nuclear reduced mass and 

y fW^ + lm y ^ (i7) 

g=u,d,s ^ q=c,b,t ^ 

for = p or n. The parameters f!j^^ are defined by 

m^4f) = {N\mgqq\N) = m,B'f\ (18) 

where [31,32] 

(jr=u,d,s 

We take the central ratios of the light quark masses from [33]: 

= 0.553, — = 18.9. 20) 
rud rud 

We take the other quark masses from [34], except for the top mass, which is taken from the 
combined CDF and DO result [35]. 

Following [36], we introduce the quantity: 



1.49, (21) 



^This expression is valid in the zero-momentum-transfer limit. For non-zero momentum exchange, the 
expression must include a form factor due to the finite size of the nucleus. See, for example, Ref. [30]. 
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which has an experimental error that is neghgible compared with others in this calculation, 
and the strange scalar density 

In terms of these, one may write 

2 + ((z-l)xy) 



(23) 



^(P) 2 X 2; - ((2; - 1) X ■ 

Proton and neutron scalar matrix elements are related by an interchange of and i?d, i.e., 

i?(") = i?iP\ Bf = B^^\ and B^^^ = Bi^\ (24) 
The TT-nucleon sigma term, S^jv, may be written as 

S.. ^\{m, + ma) x ) + sf )) , (25) 

and the coefficients /t, may be written in the forms [14]: 

- ^ - m^(l + ^) ' 

where we have dropped the (A^) superscript from /j-, and Bq. 
The value of y is related to the 7r-nucleon sigma term St^jv by 

y=l- ao/J^^r, . (29) 

The central value for ctq is estimated on the basis of octet baryon mass differences to be 
(To = 36 MeV [37-40], and the latest determination of S^^v = 64 MeV. These are the values 
assumed in our analyses of the CMSSM, NUHMl and NUHM2: the effect of varying these 
assumptions are discussed in the context of the CMSSM in [13,14]. We take this opportunity 
to reiterate the importance of measuring Stt^t as accurately as possible [14jf\ 

The spin-dependent (SD) part of the elastic x-nucleus cross section can be written in the 
zero momentum transfer limit as 

32 

asD = —GlmlA^JiJ+1) , (30) 

TT 



■^Lattice calculations arc now reaching the stage where they may also provide useful information on 
S^jv [41]. 
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where rrir is again the reduced neutrahno mass, J is the spin of the nucleus, and 



^ = -f (s('S'p) + an(5'n)) , 



(31) 



where 





V2Gf 



"27 




(32) 



The factors Aq parametrize the quark spin content of the nucleon and are only significant 
for the light (u,d,s) quarks. For definiteness, we assume 



The effects of varying these assumptions are also discussed in the context of the CMSSM 
in [14]. 



We begin with a brief discussion of detection prospects in the CMSSM. In panel (a) of 
Fig. [l| we show the {mi/2, mo) plane for tan /? = 10, Aq = 0, and > 0. The region excluded 
because the LSP is a charged f is shaded brown, and that where electroweak symmetry 
breaking cannot be obtained, resulting in unphysical < 0, in dark pink. The red dot- 
dashed contour corresponds to a Higgs mass of 114 GeV. At lower mi/2 the Higgs boson 
would be lighter, which is excluded by its non-observation at LEP [42]. We also plot a black 
dashed contour for m^± = 104 GeV, the region at lower mi/2 also being excluded by LEP. 
The green shaded region at very low mi/2 and mo is disfavored by the measured branching 
ratio for b ^ sj [43], while the light pink shaded region is favored by the measurement of 
the muon anomalous magnetic moment at the 2-a level [44]. Finally, in the turquoise shaded 
regions, the relic density of neutralinos falls within the WMAP range [45]. For the value 
tan (3 = 10 used here, the only cosmologically-preferred regions are the coannihilation strip, 
bordering the f-LSP region, and the focus-point region at large mo, where fi is small and 
the LSP is a mixed bino-Higgsino state. Over the bulk of the plane, the relic density of 
neutralinos exceeds the WMAP range by more than 2 a. There are, however, slim strips 
where ^^^"^ below the WMAP range, which lie between the strips of good relic density 
and the excluded regions they border. These portions of the plane are not forbidden, as 
there could be some additional source of cold dark matter. 

In panel (b) of Fig. [l] we show the neutralino-nucleon elastic scattering cross sections as 
functions of neutralino mass for the regions of panel (a) that are cosmologically viable (i.e., 
those where the upper limit on the relic density of neutralinos is respected), and are not 
excluded by constraints from colliders. We also plot the limits on the spin-independent cross 
section from CDMS II [46] (solid black line) and XENONIO [47] (solid red line), as well as 
the sensitivities projected for XENONIOO [48] (or a similar 100-kg liquid noble-gas detector 
such as LUX, dashed red line) and SuperCDMS at the Soudan Mine [49] (dashed black line). 
We do not show the experimental bound for the spin-dependent scattering of a neutralino 



A^f) = 0.84 ,aJ'^ = -0.43 , A^P) = -0.09. 



(33) 



3 CMSSM Models 
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tan P = 10 , p. > 




nij, (GeV) (GeV) 

Figure 1: Panels (a) and (b) show the CMSSM (mi/2,'"^o) plane and the correspond- 
ing cosmologically viable neutralino-nucleon elastic scattering cross sections as functions 
of neutralino mass for tan/3 = 10 and Aq = 0. Panel (c) shows the entire potential 
range of neutralino-nucleon cross sections as functions of neutralino mass for the CMSSM, 
with 5 < tan/3 < 55, < mi/2 < 2000 GeV, 100 GeV < ttiq < 2000 GeV, and 
— 3mi/2 ^ < 3mi/2 We consider fi < only for tan/3 < 30. Also shown are upper 
limits on the spin-independent dark matter scattering cross section from GDMS II [46] (solid 
black line) and XENONIO [47] (solid pink line), as well as the expected sensitivities for 
XENONIOO [48] (dashed pink line) and SuperGDMS at the Soudan Mine [49] (dashed black 
line). 
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on a proton as the current best limit is 0.2 pb (at m^^, = 100 GeV) [50] and is off the scale 
of our plot. The limit for the scattering on neutrons is stronger but is still only at the level 
of 0.02 pb [46]. 

When stating limits or projections, direct detection experiments assume as a prior that 
^wiMP = ^CDM, where ^cdm is that measured by WMAP. If Qwimp < ^cdm, we assume 
that the rest of the cold dark matter is due to some other source, and we rescale the plotted 
neutralino-nucleon cross sections by a factor Q^/Qcdm, so that comparison with direct 
searches is possible. 

In panel (b), we display both the spin-dependent and spin-independent (scalar) cross 
sections for neutralino-nucleon elastic scattering. In each case, there are two distinct regions 
in the (M^, a) plane, that arising from the focus-point region at < 150 GeV and relatively 
large a, and that from the coannihilation strip. In the coannihilation strip, 50 GeV < < 
400 GeV, where the lower limit on is a result of the LEP constraint on the chargino 
mass, and the upper limit on corresponds to the end-point of the coannihilation strip 
for tan/5 = 10. In contrast, the end point of the focus-point region shown is due only to the 
cut-off mo < 2 TeV that we assume. 

In addition, for mi/2 ^ 380 GeV in the coannihilation strip [m^ < 160 GeV), the nominal 
calculated mass of the lighter scalar MSSM Higgs boson is less than the LEP lower bound. 
These points are indicated by lighter shadings: lime green for the spin-independent scattering 
cross sections, and ochre for the spin- dependent scattering cross sections In general, such 
points may not be strictly forbidden, as the calculated mas^ has a theoretical uncertainty 
estimated as 1.5 GeV, and the MSSM Higgs has slightly different couplings from the SM 
Higgs for which the LEP limit was set. However, in the CMSSM the couplings are generally 
very close and (up to theoretical and experimental uncertainties) the limit should hold in 
this case. 

At very low mi/2, CDMS II and XENONIO have definitively excluded some of the region 
where rrih is below the LEP limit. We show this explicitly in panel (a), by plotting the 
reach of current and future direct detection experiments in the supersymmetric parameter 
space. Here and in subsequent parameter space scans, we display contours of the scalar 
neutralino-nucleon cross section, scaled by Q^/Qcdm if necessary, of 5 x 10^^ pb (solid green 
lines) and 10^^ pb (dashed green lines). A cross section of 5 x 10^^ pb is currently excluded 
by XENONIO for = 30 GeV and by CDMS II for = 60 GeV, and will be probed 
by SuperCDMS for up to ~ 1000 GeV. Tonne-scale liquid noble-gas detectors such as 
the proposed XENONIT or a similar detector mass for LUX/ZEP will be sensitive to scalar 
cross sections below 10~^ pb for all neutralino masses in the range 10 GeV < < a few 
TeV [48]. Indeed, they will be sensitive to cross sections below 10^^*^ pb over much of the 
preferred mass range ~ O(IOO) GeV. 

The choices tan [3 = 10 and = do not yield viable direct detection cross sections 
that are completely representative of the range of possibilities within the CMSSM. There- 
fore, in panels (c) and (d) of Fig. [l]we show CMSSM spin-independent and spin-dependent 
neutralino-nucleon cross sections, respectively, as obtained in a scan over all CMSSM pa- 

'*Note that in the focus-point region the calculated Higgs mass always exceeds the LEP lower limit. 
^We use FeynHiggs [51] to calculate the Higgs mass. 
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rameters with 5 < tan/? < 55, < mi/2 < 2000 GeV, 100 GeV < mo < 2000 GeV, 
and — 3mi/2 < ^ 3mi/2- We also allow both positive and negative /i, except for large 
tan j3 > 30, where convergence becomes difficult in the /i < case 

For future reference, we note the ranges of CMSSM cross sections for different given 
neutralino masses. At low < 300 GeV, cross sections generally exceed 10~^ pb, and 
the largest scalar cross sections, which occur for ~ 100 GeV, are already excluded by 
XENONIO and/or CDMS II [) These exclusions occur primarily in the focus-point region at 
large tan jS. On the other hand, for > 400 GeV scalar cross sections are well below 10~^ 
pb, and come from the coannihilation strip or the rapid-annihilation funnel that appears 
at large tan (3 in the CMSSM Moreover, the effective cross sections are suppressed for 
points with <^ ^cdm, and there may be cancellations at larger that suppress the 
cross sections substantially. These regions of parameter space will not be probed by direct 
detection experiments in the near future. 



4 NUHMl Models 

In NUHMl models, either /i or may be taken as a free parameter, in addition to the 
four parameters and /i sign choice in the CMSSM. To examine the parameter space, we may 
choose to fix either mi/2 or mo, along with tan/?, Aq, and the sign of /i, and perform a scan 
over the remaining two parameters [21]. 

Before exploring the NUHMl parameter space, we ffist make some general comments. 
The LSP may be either bino-like or a mixed bino-Higgsino state. If the LSP is bino-like, its 
mass is nearly equal to the bino mass. Mi, which is proportional to mi/2: m^ ~ 0.42mi/2 for 
a bino-like neutralino. Scans with fixed mi/2 can therefore provide only limited information 
on the accessibility of the full parameter space with direct dark matter detection experiments. 
When n becomes less than Mi, however, the LSP becomes increasingly Higgsino-like, and 
its mass is strongly influenced by the value of fi. Specifically, one finds that the mass drops 
off as the Higgsino component becomes more substantial. 

4.1 Exploring Large mo in Sample NUHMl Planes 

We display in Fig. |2] representative NUHMl (mA,m.o) and (/i, mo) planes in panels (a) and 
(c), and the corresponding neutralino- nucleon elastic scattering cross sections as functions 
of the neutralino mass for cosmologically viable regions of the NUHMl plane in panels (b) 
and (d), respectively. Although fixing mi/2 limits the range of neutralino masses, as already 
commented, and hence also the cross sections, the (m^,mo) and (/i,mo) planes in panels 
(a) and (c) of Fig. [2] exhibit many of the constraints and features that provide the basis 
for understanding direct detection prospects in models with non-universal Higgs masses 

^We recall that models with /i < are generally disfavoured by — 2, but we do not use this as a 
restriction on our parameter scans. 

^These results are consistent with those presented in [13]. 

^Bordering these regions, the relic density of neutralinos can be quite low with respect to flcDM such 
that some of the cross sections in panels (c) and (d) are highly scaled. 
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discussed in the following sections. In panel (a), for example, regions of the plane at very 
low rriA are excluded because the lighter scalar Higgs mass falls below the LEP bound (red 
dot-dashed line) and/or BR(6 — > sj) is too large (green shading). At large relative to 
mo, it becomes impossible to satisfy the electroweak vacuum conditions with a real value of 
fj,, leading to the pink unphysical region. The brown and black regions at very low mo are 
forbidden because they have a stau or selectron/smuon LSP, respectively. 

There are four distinct regions in panel (a) of Fig. [2] where the relic density of neutrali- 
nos falls within the WMAP range. At low m^ and mo, the selectron and stau are nearly 
degenerate, and both are lighter than the lightest neutralino, with the selectron being the 
lighter. Close to this forbidden region, there is an allowed turquoise strip where coannihila- 
tions of neutralinos with selectrons and staus bring the relic density down into the WMAP 
range. Near ttla = 420 GeV, where m^ = 2m^ as indicated by the solid blue line, a rapid- 
annihilation funnel rises out of the coannihilation strip. Between the two funnel walls where 
2m^ ~ rriA-, the relic density of neutralinos is brought below the WMAP range by an en- 
hanced annihilation rate through the direct-channel A pole. On the other side of the funnel, 
we see a continuation of the coannihilation strip, near the border of the region where the stau 
is the LSP. The final region where the relic density of neutralinos falls within the WMAP 
range occurs along the border of the unphysical region at large m^ where electroweak sym- 
metry breaking cannot be obtained, in the focus-point region. Here, as /x — ^ 0, the LSP 
becomes increasingly Higgsino-like, and annihilations to gauge bosons are enhanced. Since 
= 500 GeV in this plane, the LSP is bino-like with a mass of just over 200 GeV in all 
cosmologically viable regions except the focus-point region, where the mass of the lightest 
neutralino depends on /i and may be smaller. 

As a result, panel (b) shows cross sections possible in this plane as a strip at m^ ~ 200 
GeV, including values in both the coannihilation strip and the rapid-annihilation funnel, 
where the LSP is bino-like and m^ ~ 0.42mi/2. We note that there are points with m^ ~ 
210 GeV that have very low cross sections because they lie in the rapid-annihilation funnel 
at large mg, a possibility offered in the NUHMl because of the freedom in the Higgs masses, 
that was absent in the CMSSM. On the other hand, the focus-point region yields the largest 
cross sections. However, as the LSP becomes increasingly Higgsino-like, its mass decreases 
with fi, and the relic density eventually falls below the WMAP range. Thus, we rescale 
the cross sections at lower m^ by a factor of Q^/Qwmap, which decreases as the Higgsino 
fraction increases and Q-^ decreases. Had a larger (smaller) value of mi/2 been chosen, these 
features in panel (b) would have been shifted to larger (smaller) m^. 

The green dot-dashed line in panel (a) marks the CMSSM subspace within this NUHMl 
plane. It crosses the cosmologically viable regions in two places: at large (mA,mo) in the 
focus-point region, and at low mo in the stau coannihilation strip. The largest CMSSM 
cross sections of a few x 10~^ pb come from the focus-point region, as already seen in Fig. [T| 
whereas cross sections from CMSSM points along the coannihilation strip peak at ~ 10~^ pb. 
Moreover, they are rescaled in the portion of the CMSSM contour that extends into the area 
below the coannihilation strip, where the relic density of neutralinos can be as much as an 
order of magnitude below the WMAP range. 

Panel (c) shows an NUHMl (/i, mo) plane, again with mi/2 fixed to be 500 GeV. Many 



11 



tan P = 10 , mi/2 = 500 



> 




> 



S 1000 



niA (GeV) 



tan P = 10 , = 500 






10"^-, 




10-'- 




lo--*- 




10 = 




10'- 






ft 






10-'- 


a 






10* 




io-»- 




10-10- 




10-" 




10-'^- 





10-^-, 




10-'. 




10-*- 




10-S 




10-'- 






ft 




a 


10-'- 








10-*- 




io-». 




10-"'- 




10-"- 




10-'^- 



1000 

\i (GeV) 



spin {pass all) 
scalar {pass all) 




(GeV) 





spin (pass all) 




scalar (pass all) 




■ 1 : 


~ ' ^ - - 
1 


v*^^— 



















(GeV) 



Figure 2: Panels (a) and (c) show the NUHMl (m^,mo) and {fi,mQ) planes for mi/2 = 
500 GeV, tan/5 = 10, and Aq = (fi > in panel (a)). Panels (h) and (d) show the 
corresponding neutralino-nucleon elastic scattering cross sections for cosmologically-viable 
models as functions of neutralino mass. 
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of the features mentioned above are also present in panel (c). Note that there is a small 
region with 100 GeV< fx <200 GeV and mo < 200 GeV that may be favoured by — 2 
(light pink shading, bounded by solid black line). There are rapid-annihilation funnels and 
coannihilation strips at both positive and negative fj,. Also, at large |/i| there are unphysical 
regions where the electroweak vacuum conditions imply < 0, and the rapid-annihilation 
funnels, where rriA ~ 420 GeV, run close to these borders. There is also an analogue of 
the focus-point region in the form of crossover strips of good relic density at ~ 300 
GeV, where the LSP has a substantial Higgsino component. The BR(6 — »• 57) constraint is 
important primarily for < 0, shown in the left half-plane of (c), in which case it excludes 
the funnel and parts of the coannihilation and crossover strips. The LEP constraint on 
the Higgs mass is insignificant in panel (c), as is the upper limit on BR(i?s yU+yU^), 
which becomes important only at larger tan f3. The green dot-dashed CMSSM contours run 
roughly vertically at |/i| ~ 500-600 GeV, and in the visible part of the plane they intersect 
only the cosmologically preferred coannihilation strips at low rriQ. However, at larger mg the 
CMSSM contours would join together to form a parabola, and intersect the crossover strips 
at both positive and negative fi, outside the visible part of the plane. In the area between 
the crossover strips, the relic density of neutralinos is below the WMAP range, resulting in 
rescaled cross sections in panel (d) that decrease for the increasingly Higgsino-like LSPs at 
lower m^. There are two distinct strips of scalar cross sections in panel (d) for Higgsino-like 
LSPs: the upper one corresponding to /i > 0, and the lower to /i < 0. Cross sections for 
/i < are suppressed by cancellations in the scattering matrix element due to sign differences 
in the neutralino eigenvectors. 

We illustrate in panels (b) and (d) the reaches of current and future direct detection 
experiments in the NUHMl parameter space. As can be seen from the green dashed contours 
of asi drawn in panel (a), probing this plane will require a sensitivity below 5 x 10~* pb. 
In particular, to probe the coannihilation strip at low niA and all of the focus-point region 
would require a direct detection experiment with a sensitivity of 10~^ pb. For the {fi, mo) 
plane in panel (c), 10~^ pb would again suffice to probe the phenomenologically viable region 
between the crossover strips |^ as well as some of the coannihilation strip. Unfortunately, 
scalar neutralino-nucleon elastic scattering cross sections in the rapid-annihilation funnel 
regions of both panels (a) and (c) may be considerably lower, as these points may have very 
large values of mo - a feature not found in the CMSSM. Only a detector like LUX/ZEP with 
20 tonnes would have sufficient sensitivity to probe also the funnel. 

Fig. |3]shows corresponding (m^, mo) and (/i, mo) planes with tan P = 10 and a lower value 
of mi/2 = 300 GeV close to the best fit value for the NUHMl found in [52]. Comparing 
its panel (a) with that of Fig. [2| we see that the electroweak symmetry breaking constraint 
is more important, as is the Higgs mass constraint, whereas the 6 — > 57 constraint is less 
important. A new feature is that a region at low mo is apparently favoured by gfj_ — 2. We 
also note that the prospective reach of future searches for spin-independent dark matter 
scattering in the (m^, tuq) plane is considerably greater for mi/2 = 300 GeV than for mi/2 = 

^Within the region at very low where the chargino mass is below the LEP bound, we omit the contours 
of crsi — 10"^ pb. The relic density of neutralinos can be quite low in this region, such that 10^^ pb is not 
sufficient to probe the entirety of this already-excluded strip. 
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500 GeV, as seen by the shift in the green dashed contour between the corresponding panels 
(a) However, as seen in panel (b), there is no large increase in the spin-independent 
dark matter scattering cross section in the regions of the plane allowed by WMAP, and the 
points with larger cross sections, most of which have ~ 120 GeV or less, generally have 
rrih < 114 GeV. We also note that there are points with ~ 120 GeV and very low cross 
sections, due again to their locations at large mo along the rapid-annihilation funnel in panel 
(a). 



tan P = 10 , mi/2 = 300 




Figure 3: ^45 in Fig. but for mi/2 = 300 GeV. 

^°The dip around ttia — 280 GeV in Fig. [3|a) is due to the rescaling of the cross section in the underdense 
band between the rapid-annihilation strips. 
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Comparing the (/x, mo) plane in Fig. [Sj^c) with that in in panel (c) of Fig. [2| we again 
see increased importance for the electroweak vacuum and LEP Higgs constraints, and the 
6 — S7 constraint for /i < is also more important. We also note the appearance of a 
region at low mo that is apparently favoured hj — 2, and that the prospective reach of 
future searches for spin-independent dark matter scattering in the (mA,mo) plane is again 
considerably greater for mi/2 = 300 GeV than for mi/2 = 500 GeV. As in panel (b) of Fig. [sj 
only a few points in panel (d) escape the LEP Higgs constraint, but some of these are in the 
focus-point region with < 120 GeV, and have spin-independent scattering cross sections 
close to the present experimental upper limits. We again see points with ~ 120 GeV and 
very low cross sections, that lie along the rapid-annihilation funnel in panel (c), at large mo. 
These points have no counterparts in the CMSSM. 

After these exploratory studies introducing some properties of characteristic dark matter 
regions, we next study more systematically NUHMl planes with mi/2 and either m^i or /i as 
free parameters, keeping tuq fixed. 

4.2 The NUHMl with uia as a free parameter 

Fig. |4|^a) displays a (m^,mi/2) plane with fixed tan/5 = 10, mo = 500 GeV. The triangular 
allowed region is bounded by 6 — *• S7 at small m^, the appearance of a slepton LSP at large 
mi/2, and the absence of a consistent electroweak vacuum at larger mA and smaller mi/2- 
The diagonal blue line indicates where = mA/2: on either side there is a narrow rapid- 
annihilation funnel strip where the relic density falls within the WMAP range. This funnel 
extends only to mi/2 ~ 1200 GeV, and therefore m^ < 550 GeV for this region. There is 
another WMAP strip in the focus-point region close to the electroweak vacuum boundary, 
where the LSP is more Higgsino-like. The displayed part of the focus-point strip is cut off 
at niA = 2000 GeV, corresponding to m-^ < 600 GeV: larger values of m^ would be allowed 
if one considered larger values of itla- Apart from the region between this strip and the 
boundary, and between the funnel strips, the relic LSP density would exceed the WMAP 
range. The green dot-dashed CMSSM contour runs only through regions excluded either by 
excessive VL^h"^ or by the LEP chargino mass limit. Whilst the BR(6 — > 57) limit is also 
important at very low itla and mo, it is the constraint on the Higgs mass (shown by the red 
dot-dashed curve that is roughly horizontal at mi/2 = 400 GeV) that places a lower limit on 
the expected LSP mass for the funnel region of m^ ~ 160 GeV, where the LSP is bino-like. 

We see in panel (a) that the green dashed line indicating future sensitivity to spin- 
independent dark matter scattering is mainly in the region where the relic density exceeds 
the WMAP upper limit. SuperCDMS at Soudan would be sensitive to much of the focus- 
point region shown in panel (a), but a more advanced detector would be required if the relic 
density of neutralinos is below the WMAP range. Unfortunately, even if neutralinos make 
up all the dark matter in the universe, much of the funnel region will remain out of reach, 
even to next-generation direct detection experiments. Points associated with the funnel with 
cross sections larger than 10^^ pb fail the Higgs mass constraint. 

Panel (b) of Fig. |4] displays the scattering cross sections for the allowed points in panel 
(a). The fact that scalar cross sections in the funnel region are generally smaller than those 
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Figure 4: Panels (a) and (c) show the NUHMl (myi,mi/2) planes for ttlq = 500 GeV, 
tan/? = 10 and 50. Panels (b) and (d) show the corresponding neutralino-nucleon elastic 
scattering cross sections as functions of neutralino mass. 

in the focus-point region is reflected in the lower cutoff on the LSP mass, < 550 GeV, 
seen in panel (b) for the points with a si ^ few xlO"^" pb. In the focus-point region the LSP 
is in a mixed state rriy^ < 600 GeV, and the cross section may be much larger: by two orders 
of magnitude for the spin-independent cross section, and four orders of magnitude for the 
spin-dependent cross section. These large cross sections at large have no counterparts 
in the CMSSM, since in that model mixed states in the focus-point region have much lower 
masses. Furthermore, in the CMSSM, the focus point is reached only at large mg, whereas 
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here, mo = 500 GeV is fixed. The large scalar cross sections in the focus-point region are 
already beginning to be probed by CDMS II and XENONIO, with a few points at very low 
mi/2 already being excluded, as one can see directly in panel (b) We also note that there 
are points in panel (b) with very low cross sections even though they have < 150 GeV, 
which also have no counterparts in the CMSSM. These points associated with the funnel 
and their cross sections have been scaled down due to the low relic density in that region. 

A new feature seen in this plot is the near-vanishing of the spin-dependent cross section 
when = 350 GeV. This feature is associated with the funnel region and this group 
of points all have cross sections substantially lower than those of the focus-point region 
as remarked above. However, near = 350 GeV there is the possibility for a complete 



cancellation in Eq. (31 ) when the spin contribution from up quarks cancels that due to down 



and strange quarks for the values in Eq. ( 33 ) . While the exact position of the cancellation is 
sensitive to the values of the spin matrix elements adopted, the existence of the cancellation 
is quite robust. Of course, in this case the cross section for the scattering on neutrons will 
not exhibit a cancellation at the same place. 

Panels (c) and (d) of Fig. |4] show a corresponding analysis for tan P = 50. We first 
observe that there is no longer any substantial focus-point region, and we do not display in 
panel (c) the CMSSM contour, which lies very close to the boundary of the region where 
electroweak symmetry breaking is possible. We see that much of the plane is consistent with 
the apparent anomaly in the value of — 2. At the larger value of tan/3, the branching 
ratio for Bs /U+ /x^ is substantially larger, and the current limit [53] on the Bg — > /U+ 
branching ratio of 5.8 x 10~^ excludes points with rriA < 370 — 400 GeV (depending on the 
value of mi/2), as shown by the blue dotted curve. In this case, this limit is stronger than 
that due to the Higgs mass and excludes the lower part of the funnel and, as a consequence, 
neutralino masses < 200 GeV. 

Even though spin-independent cross sections from the funnel region are slightly boosted 
at large tan (3, much of the funnel remains out of reach of direct detection experiments, as 
this region is above the solid green contour corresponding to = 5 x 10~^ pb. Moreover, the 
absence of a focus-point region for tan jS = 50 removes a prospective source of points with 
larger cross sections. Complementary collider searches and confirmation from astrophysical 
indirect dark matter observations may be needed to discover dark matter if the relic density 
of neutralinos obtains a cosmologically-acceptable value through enhanced annihilations near 
a direct-channel pole. However, a good portion of the funnel (corresponding to acceptable 
g^^ — 2) is within the dashed green contour, and hence has a cross section above 10~^ pb. 

Fig. [5] shows how the (^^,^1/2) plane changes if we decrease tuq to 100 GeV, which 
is close to the best fit found for the NUHMl [52], from the value mo = 500 GeV studied 
in Fig. |4| The most obvious effects in panel (a) are that the region where the stau and/or 
selectron/smuon is lighter than the lightest neutralino descends from mi/2 ~ 1900 GeV 
to mi/2 ~ 500 GeV. The 6 — >■ S7 and m^j constraints also become more aggressive, with 
the result that most of the points allowed by WMAP are in a coannihilation strip where 
fni/2 ~ 420 GeV, corresponding to m^ ~ 200 GeV - which is favoured by g^ — 2, moreover. 



^^We note these points also have m^j < 114 GeV and are within the sohd green contour in panel (a) at 
very low values of mi/2 and itia- 
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There are, however, some points at both larger and smaller mi/2 in the rapid-annihilation 
funnel, and also a few points in a focus-point region close to the boundary of consistent 
electroweak symmetry breaking. We see in panel (b) of Fig. [s] that there is a wide range 
of possible values of the spin-independent cross section, particularly for small where 
nih < 114 GeV. Points with < 150 GeV and such low cross sections which are associated 
with the funnel have no analogues in the CMSSM. 




Figure 5: Panel (a) shows the NUHMl (771^,^1/2) planes for uiq = 100 GeV andtan P = 10. 
Panel (h) shows the corresponding neutralino-nucleon elastic scattering cross sections as 
functions of neutralino mass. 



4.3 The NUHMl with as a free parameter 



Fig. [6]^a) shows a mi/2) plane for tan /3 = 10 and mg = 500 GeV. We see in panel (a) that 
b excludes a region at small mi/2 for fi < 0, and that the Higgs constraint excludes 

points with low mi/2 and /x > in both the crossover strip and the rapid-annihilation funnel 
(which lies close to the electroweak vacuum boundary in this case). The intersection of the 
CMSSM contour with regions of acceptable relic density occurs only at very low mi/2, in 
regions already excluded by the Higgs constraint (and also by 6 ^ 57 for /i < 0). 

The largest cross sections excluded by BR(6 — » 57) come from points in or between the 
vee-shaped crossover strips. As the x-nucleon scattering cross sections are smaller for fi < 



than for fx > , the cross sections in these excluded regions are lower by more than an 



order of magnitude than those in regions with /i > that are still allowed. Even between the 



^■^For this reason, the upper hmit on the scalar cross section is larger for ^ > than for /x < 0, resulting 
in the lower "shadow" edge in panel (b). 
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Figure 6: Panels (a) and (c) show the NUHMl {fi, 7711/2) planes for 7nQ = 500 GeV, tan/3 = 
10 and 50. Panels (b) and (d) show the corresponding neutralino-nucleon elastic scattering 
cross sections as functions of neutralino mass. 

crossover strips, at lower where the rehc density is below the WMAP range, the scaled 
cross sections are typically a few x 10~^ pb. The largest cross sections shown in panel (b) of 
Fig. [6] come from points where the relic density is within the WMAP range (and therefore 
need no rescaling) and at low mi/2, and a few of these points are in fact already excluded 
by CDMS II. The funnel region, as has been previously noted, gives rise to rather low spin- 
independent neutralino-nucleon cross sections, even at low m^. A XENONlOO-type detector 
would be sensitive to the entire crossover strip region shown here for > 0, and much of 
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that for /i < 0, if neutralinos make up all the dark matter. We note in Fig. |6j(b) the existence 
of points with 100 GeV< <200 GeV with spin-independent cross sections ~ 10^^^ pb, 
much smaller than those found in the CMSSM for this range of masses. 

In panels (c) and (d) of Figure [6| for tan /? = 50 and uiq = 500 GeV, the stau LSP region 
has intruded so far into the parameter space that the coannihilation strip at its boundary 
joins the crossover strip to the rapid-annihilation funnel. The funnel itself is deformed 
and broadened significantly, leaving a significant region within the funnel walls where the 
relic density of neutralinos is below the WMAP range. Much of the funnel is favoured by 
(7^ — 2. The CMSSM contour crosses regions of acceptable relic density at very low /i and 
m~^, in an area excluded by 6 — S7 and the Higgs mass constraint, and again near the 
boundary of the stau LSP region, where ~ 400 GeV. The largest cross sections, which 
occur at low m^, are excluded both by 6 — 57 and by present limits on direct detection 
cross sections While some of the coannihilation strip may be probed in the future, cross 
sections in the focus-point region become as low as ~ 10"^*^ pb when the relic density is 
within the WMAP range. In panel (c), the large value of tan/? again leads to a sizeable 
contribution to Bs fi^ fi" . Points below the blue dotted curve are above the current 
experimental constraint of 5 x 10^*^ [53]. 

In the (/x,mi/2) plane for tan/3 = 10 and mo = 100 GeV, shown in Fig. ^a), a com- 
bination of the electroweak vacuum constraint, the x LSP constraint and b —>■ sj excludes 
almost all the half-plane with /i < 0. On the other hand, more points survive for /i > 0, 
in the region bounded by the x LSP constraint, the LEP Higgs constraint and the LEP 
chargino constraint. The surviving points have spin- independent cross sections below the 
current upper limits, though many points with < 400 GeV should be accessible to future 
experiments. We note that there are a few points with rrih < 114 GeV that are already ex- 
cluded by XENONIO and particularly CDMS II. On the other hand, there are some points in 
panel (b) with very low scalar cross sections which arise from regions where the relic density 
is rather low (eg. in most of the region above the WMAP strip found at mi/2 = 450 GeV) 
and the cross section has been scaled. There is also a very sharp downturn in the spin- 
dependent cross section at ~ 200 GeV with values much smaller than in the CMSSM. 
This is due to a cancellation similar to the one seen in Fig. |4] between the up, down and 
strange contributions. 



4.4 Summary 

Fig. |8] is a pair of scatter plots displaying the potential ranges of (a) the spin-independent 
and (b) the spin-dependent dark matter cross sections in the NUHMl. Comparing with 
the corresponding plots for the CMSSM in Fig. [ij we note that the spin-independent cross 
section in the NUHMl may be up to an order of magnitude larger for > 300 GeV. We also 
note the appearance of points in the NUHMl with < 200 GeV and low spin-independent 
cross sections ~ 10~^° pb. Similar features are present for the spin-dependent cross sections: 
in the NUHMl this may be ~ 10~® pb for rriy. > 500 GeV, whereas values in the CMSSM 

^■^The range of in panel (d) is restricted by the range of TO1/2 shown in panel (c). 
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tan P = 10 , m,, = 100 




Figure 7: Panel (a) shows the NUHMl (/i,mi/2) plane for ttiq = 100 GeV and tan (3 = 10. 
Panel (b) shows the corresponding neutralino-nucleon elastic scattering cross sections as 
functions of neutralino mass. 

are an order of magnitude lower. Also, the NUHMl allows the possibility of much lower 
spin-dependent cross sections for < 300 GeV than are attained in the CMSSM. 

The points with large cross sections at large are generally those with relatively large 
Higgsino components. In the CMSSM, this is possible only in focus-point regions that have 
relatively small m-^. However, in the NUHMl there are focus-point regions extending to 
larger m^, and there are also other cross-over regions where the LSP has a relatively large 
Higgsino component. Relatively low cross sections may occur in the NUHMl at points with 
relatively large values of rriQ or in the funnel regions due to the scaling applied when the relic 
density is small, as was illustrated in some of the sample planes discussed earlier. 

5 NUHM2 Models 

As already discussed, the NUHM2 has two parameters in addition to those already present 
in the CMSSM, which may be transposed into free choices of both the quantities rriA and n. 
The relatively large number renders complex a systematic survey of the NUHM2 parameter 
space. We restrict ourselves here to studies of a few parameter planes whose features we 
compare with the CMSSM and NUHMl. See [6], for further discussions of direct detection 
cross sections in the NUHM2. 
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Figure 8: Panels (a) and (b) show the entire potential ranges in the NUHMl of the scalar 
and spin- dependent neutralino-nucleon cross sections, respectively, as functions of neutralino 
mass. In both plots, we scan 5 < tan/5 < 55, < mi/2 < 2000 GeV, 100 GeV <mo < 2000 
GeV, and — 3mi/2 < < 3mi/2. The common GUT-scale value of mi = m2 is in the range 
(-2000,2000) GeV 



5.1 Sample (mi/2,^0) Plane 

We first display in Fig. [9]a sample (mi/2, mo) plane with tan P = 10 and fixed tua = 500 GeV 
and fi = 500 GeV, which reveals a couple of new features. As in the CMSSM, there is a 
region in panel (a) at large mi/2 and small mo which is forbidden because the lighter stau is 
the LSP. Just above this forbidden region, as in the CMSSM, there is a stau-coannihilation 
strip. However, jutting up from this strip at mi/2 ~ 600 GeV and ~ 950 GeV, there are 
vertical strips where the relic x density falls within the WMAP range. The double strips 
at mi/2 ~ 600 GeV form a rapid-annihilation funnel on either side of the line (indicated in 
solid blue) where m^ = 2m^. Such funnels appear only at large tan/3 in the CMSSM, but 
the freedom to choose different values of m^ in the NUHM2 permits the appearance of a 
rapid-annihilation funnel also at the low value tan/3 = 10 shown here. The other vertical 
WMAP strip appears because, as mi/2 increases relative to fi which is fixed here, the Higgsino 
fraction in the lightest neutralino x increases, which in turn increases the annihilation rate 
and hence decreases the relic density, resulting in this case in a crossover region when mi/2 ~ 
950 GeV. At slightly higher mi/2 = 1020 GeV, there is a rapid coannihilation pole (through 
Z-exchange) between the lightest and next-to-lightest neutralinos resulting in a narrow region 
with suppressed relic density and hence a suppressed scalar cross section. 

These novel regions are clearly visible in panel (b) of Fig. [9] The elastic scattering cross 
sections do not vary rapidly as m^ increases, until the funnel region at m^ ~ 250 GeV is 
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Figure 9: Panel (a) shows the NUHM2 (mi/2,mo) plane for tjia = 500 GeV, /i = 500 GeV 
and tan/3 = 10. Panel (b) shows the corresponding neutralino-nucleon elastic scattering 
cross sections as functions of neutralino mass. 



reached. The vee-shaped suppression in the cross section arises from the increasing value 
of mo as one rises up the funnel. The two sides of the funnel approach each other as mg 
increases, eventually joining together and resulting in a minimum value ctsi ~ 5 x 10^^^ pb 
where the two sides of the funnel meet (at a value of mo > 2000 GeV, and hence invisible in 
panel (a) of Fig. [9]). After this excitement, the cross section continues to rise gradually as 
one follows the coannihilation strip, until the crossover strip is reached at m^ ~ 400 GeV. 
Here the cross section decreases again as mg increases, to values even smaller than in the 
rapid-annihilation funnel, before rising again and finally declining towards the end of the 
coannihilation strip. Note that the entire region to the right of the transition strip is viable, 
albeit with a relic density below the WMAP range. The cross section will thus be reduced 
due to scaling. Because the neutralino is predominantly a Higgsino here, its mass is given 
by /X rather than mi/2 and so points at large end in panel (b) because of our choice of 
fixed ji. 



5.2 Sample (m^,mi/2) Planes 



Fig. 10 shows two NUHM2 (m^,mi/2) planes. We see in panel (a) for mg = 500 GeV, 
/i = 200 GeV and tan /5 = 10 a strip with the relic density in the WMAP range that extends 
to large mA at mi/2 ~ 290 GeV. This is the transition strip, and above it, the relic density is 
always below the WMAP range. There are other strips below this, but they fall in a region 
where the chargino mass is below 104 GeV. As seen from the location of the dash-dotted red 
line, all of these cosmologically-preferred strips have mh < 114 GeV. In this case, the funnel 
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region shown by the sohd blue hne where tua = 2m^, occurs past the transition region and 
there rapid annihilation further suppresses the relic density. It is clear from the shape of 
this line that ~ /i = 200 GeV for larger mi/2- We see in panel (b) that the low-m/j 
points typically have spin-independent cross sections ~ 10~^ pb, and are largely excluded 
by the XENONIO and CDMS II experiments. On the other hand, essentially all the points 
with rrih > 114 GeV survive the direct dark matter search experiments, so far, in particular 
because the effective cross section is suppressed by the small relic density. 




Figure 10: Panels (a) and (c) show the NUHM2 (^^,^^11/2) planes for rriQ = 500 GeV, 
fi = 200 GeV, tan/? = 10 and 50. Panels (h) and (d) show the corresponding neutralino- 
nucleon elastic scattering cross sections as functions of neutralino mass. 
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Turning to the corresponding (myi,mi/2) plane for tan/3 = 50, we see a large region at 
large tjia and/or mi/2 that is excluded because the LSP is charged. There is also a large 
region at smaller mi/2 that is excluded by the 6 — > S7 constraint . At large tan/?, the 
constraint from Bg f^'^ fJ-' is important, excluding values of ^ 300 GeV in this case, 
as shown by the dotted blue curve. Once again, although the funnel region, shown by the 
solid blue line where tjia = Sm,., has rrih > 114 GeV in between the shaded regions and 
to the right of the —>■ ^'^ ^~ constraint, the relic density is small as this occurs past the 
transition region (which here is excluded by 6 ^ 57) and the neutralino has a large Higgsino 



component. As seen in panel (d) of Fig. 10, all the allowed points have ^11 = 200 GeV. 
Typical spin-independent cross sections are somewhat larger than in the case tan/3 = 10, 
shown previously in panel (b). 



We see in Fig. 11 the evolution of these features when /i = 500 GeV, with the other 
parameters left unchanged. For tan/? = 10 in panel (a), we see that the transition strip at 
mi/2 ~ 950 GeV is now split by the two halves of a rapid- annihilation funnel. In contrast to 



Fig. [10} a portion of the funnel now lies below the transition strip and the WMAP density 
can be realized. Above the crossover strips, at larger mi/2, the LSP becomes Higgsino-like. 
Only small portions of the WMAP strips have uih < 114 GeV. We see in panel (b) that 
the scatter plot of the spin-dependent and spin-independent cross sections has a feature at 
~ 430 GeV, corresponding to the opening out of the rapid-annihilation funnel. The 
points at larger have relatively large spin-independent cross sections, higher than in the 
CMSSM for similar values of m^, reflecting the Higgsino nature of the LSP in these cases. 



Turning now to panel (c) of Fig. 11 , for tan /? = 50, we see again the increased importance 
of the neutralino LSP and 6 ^ 57 constraints and the constraint from Bg — > /i"*" /i~ excluding 
small values of itia- We also see the WMAP strips on either side of the solid blue line where 
niA = 2?7i^, and notice that although almost all of them have m/j > 114 GeV only the upper 
part of the funnel is allowed by Bg fi^ . The would-be transition strip has been pushed 

into the region with a charged LSP. The corresponding values of ~ 250 GeV, as seen in 



panel (c) of Fig. 11, and the spin-independent cross section may be as large as ~ 10~^ pb, 
namely considerably larger than in the CMSSM, but decreasing for points with a low relic 
density. 

5.3 Sample (/i,mi/2) Planes 



Fig. ga) displays the NUHM2 (/i, mi/2) plane for mo = 500 GeV, niA = 500 GeV and 
tan/? = 10. Highly visible at large mi/2 are regions excluded by the neutralino LSP require- 
ment. In the dark (very dark) shaded regions, the LSP is the stau (selectron/smuon). Much 
of the remaining fi < half-plane is excluded hj b —>■ s'j, while the lower part of the /i > 
half-plane with mi/2 < 300 GeV has mh < 114 GeV. Apart from a very small section of 
a crossover strip near fi ~ —500 GeV and mi/2 ~ 1200 GeV, the WMAP density range is 
attained only in the /i > half-plane. This occurs in a rapid-annihilation funnel on either 
side of the blue line where tua = 500 GeV = 2m^, and along a crossover strip extending to 
higher (and somewhat lower) mi/2- 



Panel (b) of Fig. 12 displays the corresponding dark matter scattering cross sections. We 
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Figure 11: Panels (a) and (c) show the NUHM2 {1^^,1711/2) planes for itlq = 500 GeV, 
/i = 500 GeV, tan (3 = 10 and 50. Panels (b) and (d) show the corresponding neutralino- 
nucleon elastic scattering cross sections as functions of neutralino mass. 

see that the spin-independent cross section may be as high as 5 x 10~^ pb for ~ 500 GeV, 
to be compared with a maximum of ~ 10"^ pb in the CMSSM: this possibihty occurs for 
points towards the top of the crossover strip. At small m^, there are a few points with 
rrih < 114 GeV that are excluded by the XENON 10 and CDMS II experiments, which 
occur at the bottom end of the crossover strip. In between, for ~ 250 GeV, there is 
a suppression of the maximum spin-independent cross section, corresponding to the rapid- 
annihilation funnel extending to large /i. Foreseen experiments should be able to cover all 
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Figure 12: Panel (a) shows the NUHM2 (/x,mi/2) plane for uiq = 500 GeV, = 500 GeV 
and tan/3 = 10. Panel (b) shows the corresponding neutralino-nucleon elastic scattering 
cross sections as functions of neutralino mass. 



the WMAP strip except a small portion of this funnel. On the other hand, there are many 
points with lower effective cross sections, suppressed by the low relic density. These include 
some at low < 200 GeV with lower cross sections than those found in the CMSSM. 

Fig. [I3| displays similar NUHM2 (/i, mi/2) planes for mo = 1500 GeV, = 1000 GeV and 
tan (3 = 10. In this case, we see in panel (a) that the neutralino LSP and b ^ s'j constraints 
have no effect. Almost all the half-planes for both signs of fi have mt > 114 GeV, and the 
rapid-annihilation funnel has risen to mi/2 ~ 1100 GeV, corresponding to the higher value 
rriA = 1000 GeV. In addition to the rapid-annihilation funnels, almost all the crossover strips 
are allowed for both signs of /x. Panel (b) has features rather similar to those in panel (b) 
of Fig. 12, with relatively large spin-independent cross sections ~ 10~^ pb possible for all 

< 800 GeV, along the crossover strips, and dips around = 500 GeV, corresponding 
to the rapid-annihilation funnels. As in the previous figure, low effective cross sections are 
again possible for points with low m^, if they have a suppressed relic density, as occurs 
between the two crossover strips for the different signs of fi. 



5.4 Summary 



Fig. 14 displays the potential ranges of (a) the spin-independent and (b) the spin-dependent 
dark matter scattering rates in the NUHM2. Comparing with the corresponding plots for 
the CMSSM in Fig. [H we note that the spin-independent cross section in the NUHMl may 
be up to an order of magnitude larger for > 300 GeV. As discussed in the previous 
subsections, the neutrahno LSPs at these points typically have large Higgsino components. 
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Figure 13: Panel (a) shows the NUHM2 (/i, mi/2) planes for ttiq = 1500 GeV, ttia = 1000 
GeV and tan f3 = 10. Panel (b) shows the corresponding neutralino-nucleon elastic scattering 
cross sections as functions of neutralino mass. 

despite their large masses, and lie along crossover strips. This feature is common with the 
corresponding scatter plots for the NUHMl shown in Fig. [8] We also note the appearance 
of NUHM2 points with < 200 GeV and low spin-independent cross sections ~ 10"^'' pb. 
These points are typically in regions between the crossover strips, in regions with a very 
suppressed relic density. Similar features are present for the spin-dependent cross sections: 
in the NUHM2 this may even be ~ 10~^ pb for > 500 GeV, whereas values in the 
CMSSM are over an order of magnitude lower. Also, the NUHM2 allows the possibility of 
much lower spin-dependent effective cross sections for < 300 GeV than are attained in 
the CMSSM. These features are again similar to those found in the NUHMl. 



6 Conclusions 

We have explored in this paper the possible ranges of dark matter scattering rates in the 
NUHMl and NUHM2, and compared them with the ranges attainable in the CMSSM. 
We have seen that the ranges that could be found in the NUHMl and NUHM2 are both 
significantly broader than is possible in the CMSSM. In addition to the expectations for 
direct detection from broad scans, we have displayed calculated cross sections in specific 
slices of the NUHMl, 2 parameter spaces to highlight the physical processes in x — p elastic 
scattering. Larger cross sections are possible at large m^, and smaller cross sections are 
possible at small m^. The high-mass points with the largest cross sections occur typically 
for points in the crossover regions with relatively large Higgsino components. Mixed Bino- 
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Figure 14: Panels (a) and (b) show the entire potential range of neutralino-nucleon cross 
sections as functions of neutralino mass for the NUHM2 for the scalar and spin dependent 
cross sections respectively. In both scans, we scan 5 < tan/? < 55, < mi/2 < 2000 GeV, 
100 GeV < mo < 2000 GeV, — 3mi/2 < < 3mi/2, and the GUT-scale values of m\ and 
m2 are each in the range (—2000,2000) GeV. 



Higgsino states with the a rehc density in the WMAP range may occur at much larger masses 
than those found in the focus-point region of the CMSSM. Several mechanisms yield points 
with small effective cross sections. Some points have a very suppressed relic density, e.g., 
points with even larger Higgsino components that are located between crossover strips. In 
other cases, the cross section may be suppressed because the point is in a rapid-annihilation 
funnel region at larger mo and smaller mi/2 than is possible in the CMSSM, or in a crossover 
region at larger mo, or at larger m^ or fi. 

Present direct dark matter searches for spin-independent scattering, in particular the 
XENONIO and CDMS II experiments, are beginning to chip away at the parameter spaces 
of the CMSSM, NUHMl and NUHM2. However, most of the points excluded so far have 
small m^ and m^ < 1 14 GeV. On the other hand, prospective searches for spin-independent 
scattering could be sensitive to NUHMl or NUHM2 models with LSP masses as large as 
800 GeV, the largest we have sampled, whereas they would be sensitive only to m^ < 350 GeV 



within the CMSSM Thus, observation of spin-independent scattering by a heavy LSP 
would be a good diagnostic for a more complicated supersymmetric model than the CMSSM, 
and perhaps a hint for a neutralino with a significant Higgsino component in one of the 
crossover regions. Conversely, if some other experiment, e.g., at the LHC, establishes the 



^'^Unfortunately, searches for spin-dependent dark matter scattering are still some distance away achieving 
the sensitivity needed to probe the classes of models discussed in this paper. 
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existence of a light LSP, but its spin-independent scattering is not seen, this could also be 
a hint for some model more complicated than the CMSSM in which the relic density is 
suppressed, e.g., by one of the mechanisms described in the previous paragraph that operate 
in the NUHMl and NUHM2. 

Neither the LHC experiments nor direct dark matter searches can, by themselves, tell us 
all we would like to know about the manner in which supersymmetry is broken. However, 
as we have illustrated in this paper, the combination of LHC experiments and direct dark 
matter searches could provide some interesting hints. 
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